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SUMMARY 

Approximate solutions are obtained describing mainstream flows con- 
fined to regions of small angle variation over flat surfaces for three- 
dimensional , laminar, incompressible, thin boundary- layer flows having 
similarity with respect to stationary polar coordinate systems. The solu- 
tions, summarized in a table, include accelerating or decelerating flows 
and stagnat ion- point , spiral, or circular flows. An experimental compari- 
son of limiting overturning at the wall showed good agreement for the 
first 10° of turning of circular mainstream flow. 


INTRODUCTION 

In addition to providing an insight into secondary- flow behavior 
associated with laminar boundary- layer flows, the experimental investiga- 
tions of references 1 to 3 demonstrate that the information thus obtained 
for laminar flows can be used to interpret and to correlate flow measure- 
ments taken in turbomachines at operational conditions. These experimental 
investigations thereby provide an important link between applied turbo- 
machine research and the similarity- type boundary- layer analyses developed 
in references 4 to 12. The link is further strengthened by the combined 
theoretical and experimental investigation of reference 13. In reference 
13, boundary- layer similarity solutions are obtained for main flows con- 
sisting of streamline translates (i.e., the entire streamline pattern can 
be obtained by translating any particular streamline parallel to the 
leading edge), and the theoretical predictions of boundary- layer over- 
turning (more than mainstream turning) near the surface are in close 
agreement with experimental results obtained by tracing the boundary- 
layer streamlines with smoke flow-visualization techniques. 
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Using a generalized similarity variable tj, reference 14 extends 
these results analytically to obtain all possible flows with boundary 
layers having classical similarity with respect to stationary rectangular 
coordinates. The dimensionless boundary- layer velocity components in the 
plane of the surface are assumed to have similarity with respect to their 
respective coordinates. This similarity is expressed by means of two 
suitably defined functions of the similarity variable. The boundary- layer 
equations are then transformed to equations involving the mainstream flow 
components, their derivatives, the similarity functions, and their deriva- 
tives. All the mainstream flows are then determined for which the trans- 
formed boundary- layer equations reduce to ordinary differential equations 
in the similarity functions and their derivatives. Four distinct fami- 
lies of such mainstream flows are obtained in reference 14, including 
cases of accelerating or decelerating flows for quite general streamline 
paths. The main-flow streamlines are not required to be translates, nor 
are they restricted to regions of small turning in reference 14. 

In reference 15, solutions are obtained for the mainstream flows 
whose boundary- layer velocity components in the plane of the surface have 
similarity with respect to the corresponding polar coordinates. Thus, 
exact solutions are obtained for spiral, circular, and stagnation- point 
flow configurations with no restrictions on mainstream turning. For the 
solutions thus obtained, however, a proper leading edge cannot be defined. 
(A proper leading edge, which corresponds theoretically to a real physical 
leading edge, would be a line or curve of zero boundary- layer thickness 
on the surface, downstream of which the boundary layer develops.) 

The present investigation extends the analysis of reference 15 by 
considering the flows in a sector-region of small central angle 0. The 
purpose of this investigation is to determine mainstream flow solutions 
for which the transformed boundary- layer equations reduce to ordinary 
differential equations. Solutions are obtained for four new families of 
mainstream flows with boundary layers having similarity with respect to 
the polar coordinates in the plane of the surface. Included here are 
cases for flows over well-defined leading edges. It is important to note 

that experimental investigations (refs. 1 to 3 and 13) indicate that in 
typical turbomachine configurations a large portion of the end-wall 
boundary layer at the inlet to a passage has completely crossed from the 
pressure to the suction side of the passage when the mainstream has been 
turned less than 30°. Thus, in describing physical flow, it appears not 
unreasonable to restrict the analysis to small central-angle sectors. 

The regions where this assumption might be considered reasonable are 
established by a theoretical and experimental comparison of the boundary- 
layer limiting flow deflection (ref. 13) in a circular two-dimensional 
channel* 
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SYMBOLS 

a,b,c,C 

constants 

F,F(t]) 

function of similarity parameter, u = UF ! (r)) 

«-T]) 

function of similarity parameter, eq. (45) 

f(r),f*(r) 

arbitrary functions of r 

/(r,e) 

a.:rb ±"fc:rajry function of 1 t and Q 

G,G(tj) 

function of similarity parameter, w e WG’(ti) for W ^ 0, 
w s WG* (ti) for W = 0 

g,g(r,0) 

function of coordinates r and 9 

h(0) 

function of 9 

k,m,n 

constants 


cylindrical coordinates 

t 

constant 

u,w 

mainstream velocity components in 6 and r direc- 
tions, respectively 

u,v,v 

boundary- layer velocity components in 9, y, r direc- 
tions, respectively 

W,W(r,0) 

function of coordinates r and 9, w = WG'(ti) for 
W = 0 

r 

boundary- layer deflection at surface 

n 

similarity variable t] = yg(r ,9)/\jv 

V 

coefficient of kinematic viscosity 


Subscripts: 


i = 1,2,3, . . . 

index numbers 

Superscripts: 



Primes denote differentiation 
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ANALYSIS 

Boundary- Layer Equations in Stationary Cylindrical Coordinates 

The three-dimensional laminar incompressible thin boundary- layer 
equations in cylindrical coordinate form for flows over flat (or nearly 
flat) surfaces with stationary coordinate axes as shown here 



are given by 


uw u du Su 5u S^u UW U SU TT SU 

T 4 -735 + ' , Sf + v 37- v ^ = T + 73e + w 5F 
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in the tangential direction and 


n 

XT U Sw &W 

r + r ^0 + W "5r + V cty 




r 


U dw 
+ r^e 


+ w 


dW 

3? 


(lb) 


in the radial direction, where u, w, and v are the boundary- layer 
velocity components in the 0, r, and y directions, respectively. Con- 
sistent with the restriction to thin boundary- layer flows over flat (or 
nearly flat) surfaces as required for the formulation of the boundary- 
layer equations (eq. (l)), the mainstream velocity components are 


U = U(r ,e) 

(2a) 

W = w(r,e) 

(2b) 


The equation of continuity for the boundary- layer flow is 

1 du 5w w bv ~ 

rSe + 3? + 7 + ^ =0 

The appropriate boundary conditions are 

u = w = v = 0 for y = 0 
u -+ 

> as y 
w W I 


(lc) 


(l d) 

(le) 


Similarity with Respect to Stationary Polar Coordinates 

The boundary- layer equations may be transformed by the use of a 
generalized space variable 



Tt «-*= g(r,0) 

V v 

(3) 

and by defining 

u ■ UF'(n) 

(4a) 


w = WG 1 (t] ) 

(4b) 

to a new system of 

coordinates r, 9, and t). 

The definitions (4a) and 


(4b) are the requirements for similarity of the boundary- layer velocity 
components in the plane of the surface with respect to their corresponding 
polar coordinates . 
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The choice of polar coordinates (as in ref. 15) requires an additional 
precaution beyond those needed for similarity with respoct to rectangular 
coordinates. In rectangular coordinates when either U or W is zero 
(see ref. 14), the mainstream flows are straight, there is no secondary- 
flow overturning in the boundary layer, and complete similarity solutions 
have been obtained for the equations of the resulting two-dimensional 
boundary- layer flows (refs. 4 and 5). In the present case, however, when 
the mainstream radial component W = 0, there is curvature of the main- 
stream flow, U ^ 0, and three-dimensional boundary- layer overturning 
results (i.e., w ^ 0, for W = 0). Under these conditions, equation (4b) 
does not apply. Instead, a new function W(r,0) is defined for 

w — WG'(ti), W ^ 0 (4c) 

Accordingly, it is convenient to treat W / 0 flows separately from 
W = 0, W ^ 0 flows. 

W z/: 0. - When the mainstream flow has both U and W components, 
the corresponding boundary- layer velocity components are defined by 
equations (4a) and (4b) as functions of the similarity parameter T). 
Corresponding "to "the conditions of* no flow at the surface (aq» (id)), the 
boundary conditions on F 1 and G ! are 

F 1 (0) = G f (0) = 0 (5a) 


Corresponding to the condition (le) that the u and w boundary- layer 
components merge smoothly into main- flow components U and W, 
respectively, 

lim F 1 (t| ) = 1 

T) 00 

lim G* (t] ) = 1 

TJ “► °° 

Now v may be determined by integration of the continuity equation 
using (4a) and (4b): 



v = 


g 


& 


dU 

30 


U d 


In 

3T 


! ) p+ ( 


Sw 

"5r 


W 


5 In g 
5r 



U a In g F , 
r d0 y 


- W 


a In g 

ar 


yG' 


+ /(r,0) 


where /(r ,0) is an arbitrary function arising from integration. 


( 6 ) 
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In order that v = 0 for y = 0 as required, it is possible with- 
out loss of generality to set the boundary conditions 

F(0) = G(0) = 0, and /(r,0) = 0 (5c) 

(See appendix C, ref. 14, for a discussion of the necessary and sufficient 
boundary conditions. ) 

Upon substitution of equations (4a), (4b), and (6), equation (la) 
becomes (for W ^ 0) 

| (F'G 1 - GF" - 1) + ^ (F< 2 - FF” - 1) + W (F»G* - l) + 

JL d l n g 2 FF » + W d ^2. GF" - GF" - g 2 F'" =0 (7) 

2r 50 2 dr dr 

and equation (lb) becomes 

£ (X - F' 2 ) + 2 (F'G' - i) ♦ g (O' 2 - GO" - 1) - 7 £ *r * 

U d I q-.e!- + | Lg-s! GG" - | GG" - g 2 G’" = 0 (8) 

As in references 14 and 15, the purpose of this investigation is to 
determine mainstream flow solutions for which the transformed equations 
(7) and (8) reduce to ordinary differential equations. As an extension 
of reference 15,* the present analysis considers flows restricted to a 
sector-region of small central angle 0. The mainstream flow conditions 
are sought which make the coefficients of the functions of t] 
proportional. 

The most general approach would be to rewrite (7) and (8), grouping 
the coefficients of like terms in G,F and their derivatives, and then 
to require proportionality of these grouped coefficients. It can be 
shown, however, that no cases arise beyond those obtained more simply by 
requiring proportionality of the individual coefficients in (7) and (8). 
Under these ordinary differential equation conditions (abbreviated to 
o.d.e. conditions), the common variable terms in the equations may be 
divided out, leaving ordinary differential equations for F and G. The 
actual numerical solutions of the ordinary differential equations are not 
attempted herein. 
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For convenience, the coefficients for the functions of T) in equa- 
tions (7) and (8) are presented here in the order of their appearance. 
With W ^ 0, they are 


CD 

w 

r 


© 

dw 

© 

1 du 

r 

© 

g 2 

© 

w d 

In U 
dr 


Wr 

© 

u d 
2r 

In g 2 

~5e — 


U d In W 
r W~ 

© 

W d 
2 

In g 2 

3? 




The o.d.e. conditions require these nine coefficients to be pro- 
portional to each other. 


W = 0 . - When W = 0, the corresponding boundary- layer equations 
(la) and lb ) become 


uw u du du du o u U dU 

r + r ^8 + w 3r + V Sy V ^y2 r ^0 

u 2 u dw dw dw d 2 w U 2 

r + r^0 + W ^r + V "Sy V ^ y 2 “ r 


(9a) 


(9b) 


The equation of continuity for the boundary- layer flow remains unchanged: 

1 du dw w dv 


r^ + ^ + r + 3y - ^ 

The boundary conditions now are 

u = w = v = 0 for y = 0 
u U 

w -> 0 


(lc) 


as y 


(9c) 




4234 


NACA TN 3890 


9 


For main flows such that W = 0, define 

u s UF* (i]) 
w s WG ’ ( i] ) 


(4a) 

(4c) 


where W = W(r ,Q) ^ 0. The boundary conditions on F* and G l required 
to satisfy boundary conditions on u and w in equations (9) (W = 0, 

W(r ,9) ^ 0) are 


F'(0) = G‘(0) = 0 
lim F 1 ( rj ) = 1 

T| -*■ oo 

lim G'(ti) = 0 


(10a) 

(10b) 


The expression for v obtained by integration of the continuity 
equation (lc) is the same as (6), with W being replaced by W: 


v = 


- Vv 


6 


dU U d In g 




"5F 


) F + (I 


_ y d In g 


Sr 




U d In g — d In g f, 

7 -§q yP' - W — g - r & yG* + /(r,6) 


( 11 ) 


As before, the boundary conditions chosen as sufficient to provide 
that v = 0 for y = 0 are 


F(0) = G(0) = 0 
/(r,0) = 0 

Substitution of equations (4) and (ll) into (9a) and (9b) produces 
^ (F‘G* - GF") + ^ || (F‘ 2 - FF" - 1) + W U F*G* + 


(10c) 


U d In e 2 


W 5 In g 


~5F 


dw 


£_ 

Wr 


_ - 5 ^-=- FF" + f -jj-*- GF" - £ GF" - g‘F"' - 0 
(1 - F' 2 ) + ? P'G* + f (G - 2 - GO”) - ± g W + 


( 12 ) 


U d i n g 2 PG" + | GG" - W GG" - g 2 G "' = 

2r 60 2 dr r 


0 


(13) 
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The argument concerning determination of the o.d.e. conditions by 
means of relations between the coefficients of the functions of t] in 
(12) and (13) remains unchanged. These coefficients, it may be noted, are 
the same as the coefficients for equations (7) and (8), respectively, with 
W replaced by W. 


As the flow is considered restricted to regions of small variation 
of 0 (i.e., regions that are narrow sectors having small central angles), 
coefficients will be neglected that are of second and higher order in 0 
relative to the other coefficients. In order to do this, it is assumed 
that U(r,0), w(r,0), and W(r ,0) are expressible as 


and that, in the region of interest chosen for convenience about 8=0, 
^(0) and h 2 (0) are adequately represented by 


Hence, U, W, and W in the neighborhood of 0=0 are considered to be 
defined by 


As will be seen later, in some cases this assumption further results in 
solutions that have properly defined leading edges. Solutions are ob- 
tained for the main flows in the following manner. 


Solutions for Small Variation of 0 


U(r,0) = fj(r)h 1 (0) 

W(r , 0) = f*(r)h2(0), W / 0 
W(r,0) = f|(r)h 2 (0), W = 0 


h x (0) * a0 t 
h 2 (0) * b0 m 



W(r,0) = f 2 (r)0 m J 
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By substitution, the coefficients of the functions of F and G 
in equations (7) and (8) for W ^ 0 or in equations (12) and (13) for 
W = 0 become 


fo (r ) 

© 8" 




© | f 1 (r)0 t - 1 
fp(r ) 

© f^7 

q f d .^s 2 e t 


© 

© 

© 

© 


2r 

80 

f 2( r ) a 

In g 2 

2 

dr 

f£(r)e m 


s 2 


(iiW) ! 

5 

- e 2t_ 


e m } 


m 


© m - 0 t_1 


J 


(15) 


As before (refs. 14 and 15), the objective is to find conditions 
(o.d.e. conditions) that make these coefficients proportional to one 
another. Then the common r and 0 factors may be divided out, and the 
transformed equations reduce to ordinary differential equations. Here, 
however, with 9 small, it is assumed that terms of second and higher 
order in 6 relative to the rest are negligible. 

The procedure will be to establish the relations between constants 
m and t (the powers of 0 involved) or the assumptions concerning 
f-^(r) or f2( r ) that will lead to proportionality among the terms. The 

following possibilities are suggested. 
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m = t - 1 . - When m = t - 1, 8 m 
coefficients listed in (15) except (7) 


is a common factor of all the 
and may be divided out, yielding 


©T 

m + 1 


© 


F~ f i 


f 2 f i 






f l S In g 2 


2r 

£2 d In g 2 
2 dr 


© q 

© g 2 0' m 
*2 


© 


1 e 2 


rf. 




m 


J 


(16) 


Coefficient © is the coefficient of the F'" and G"' terms in the 
transformed boundary- layer equations. It must not he permitted to vanish, 
for that would reduce the order of the transformed equations, and the 
number of boundary conditions on F and G would then exceed the order 
of the equations. Consequently, © is made proportional to ©, with 
the result that 

g 2 = a^‘0 111 (17) 

Coefficient © is therefore independent of 9. Coefficient © is of 
second order in 9 relative to the rest and is therefore neglected. 

From o.d.e. conditions on © and ©, 


f l ~ a 2 f 2 


(18) 
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and, from © and ©, 


f„ = a_r 


(19) 


Therefore, combining ( 14 ) , (17), (18), and (19) and redefining the 


constants a^ for convenience results in 


U = ar n e m+1 
W = br n 0 m , W / 0 




g 2 = cr n-l 0 m 


( 20 ) 


J 


The corresponding ordinary differential equations are obtained by 
substitution in equations (7) and (8): 


b(n + l)(F*G l -l)- b ^ 2 + ^ GF" + a(m + 1) [(F« ) 2 - l] - FF" - cF"' = 0 

( 21 ) 

am(F'G« - l) + bn [(G* ) 2 - l] - GG" - pQ" . C G’" = 0 

( 22 ) 


with equation (5) giving the boundary conditions. When W = 0, 

W = br n 0 m ^ 0; and substitution into equations (12) and (13) yields 

b(n+ l) F*G* - b GF" + a(m + l) JjF 1 ) 2 - lj - a ^ m t _) FF" - cF'" = 0 

(23) 

am(F*G‘) + bn(G') 2 - GG" - pg" - cG"' = 0 (24) 

with equation (10) giving the boundary conditions. It is important to 
note that the boundary conditions on G' are different for W ^ 0 in 
equations (2l) and (22) from those for W = 0, W ^ 0 in equations (23) 
and (24). While a more general case appears to be 


m = t - k, k = 1,2,3 
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only the case 
this analysis 
compared with 

m = t + ; 
all terms by 


From the 
has the form 


discussed here, k = 1, has any practical significance for 
When k ^2, then, for small 0, U is negligible when 
W, and the resultant flows are cases of U = 0. 

L. - When m = t + 1, substituting into (15) and dividing 

0 m -2 gj_ ve 


0 -# qZ 




m - 1 4, 

r f l 


0 | qe 2 

, 0 . f l d In g 2 

© 27^0 — e 


©^ 


f 2 d In g 2 .2 


2 5r" 


f'Q c 


0 g 2 0 2_m 


© 


(fi ) 2 


rfc 


B fl 

r J- 


J 


(25) 


o.d.e. conditions on (j) and (§) it can be seen than g c 


„Z i j? /^in - 2 

g =~ f ± e 


( 26 ) 
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Coefficient (4) is therefore independent of 6 and so ©, ©, (5) , and 
are considered negligibly small. The modified coefficients are now 
written 


© 0 


(D 



0 o 

© m - 2 ^ 

2 r~ f l 

© 0 > 
© 0 



© 



J 


(27) 


From the o.d.e. conditions © and © it can he seen that 


f l = b 2 f 2 


(28) 


and no further restrictions on the form of f^ or fg are required, as 
all coefficients are proportional. Then using (14), (26), and (28), U, 
W, and are written 


U = af(r)e m_1 

"N 



W = bf(r)0 m , W / 

0 

> 

(29) 

g 2 = 7 f(r)e m_2 




r 

J 



and equations (7) and ( 8 ) become 




a(m- 1) HfM 2 - l] -f 

FF 

" _ cF"' = 0 

(30) 
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- £( F ') 2 - i] + am(F'G' - l) - ™ PG" - cG m = 0 (3l) 

The boundary conditions are given by equation (5). 

When W = 0, W = bf(r)0 m ^ 0, substitution into equations (12) and 
(13) produces 

a(m - 1) jjF‘) 2 - lj - ™ FF" - cF m = 0 (32) 

- [(F‘) 2 - l] + am(F , G‘) - ™ PG" - cG'" = 0 (33) 

Again, while 


m = t + k, k = 1,2,3 

appears more general, only k = 1 actually applies here. When k ^ 2, 

W = 0 under the present assumptions, and the resulting flow is equiva- 
lent to taking b = 0 in equations (20) or (29). 

m = 0 . - When m = 0, from o.d.e. conditions (l) and (8) and the as- 
sumption that terms of second or higher order in 0 relative to other 
terms may be neglected, t = 1, -1, or 0. The case t = 1 corresponds 
to the case of equation (20) and t = -1 to the case of equation (29). 

If t = 0, U = U(r), fj is proportional to f^, and therefore 

U = ar n 
W (or W) = 


This case and the resulting ordinary differential equations are the same 
as were obtained by taking U = U(r), m=0 in the exact solutions 
(ref. 15). 

t = 0 . - When t = 0, examination of the o.d.e. conditions discloses 
no cases not already obtained by the analysis. 

f^(r) = a or f2( r )' = h. - Where fp(r) is taken to be a constant, 
a, U = U(0) alone and the coefficients in (15) become 
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© 

© 

© 

© 

© 

© 


t 0 t-l 
r 

0 

1 d In g 2 n t 

27 7e 0 

f 2 5 In g 2 

T — 0 

f«e m 

6 2 

1 0 2t-m 




© 


m 0 t-l 
r 


J 


(35) 


Just as before, two significant possibilities occur, m = t - 1 and 
m = t + 1. When m = t - 1, equation (18) applies and f 2 is a constant. 

This corresponds, of course, to n = 0 (eq. (20)). When m = t + 1, the 
analysis leading to equation (28) applies, and the resultant flows cor- 
respond to f(r) constant in equation (29). 

When f 2 (r ) is constant, the same results are obtained and f^(r) 
must likewise be constant. 


RESULTS AND DISCUSSION 

The analysis of three-dimensional, laminar, incompressible boundary- 
layer flows having similarity with respect to polar coordinates has led 
to two different categories of solutions. In the first are the solutions 
for mainstream flows described in reference 15. In the second category 
are perturbation- type solutions obtained here whose validity is restricted 
to regions of small variation of angle 0. The mainstream flows for the 
latter category are described by equations (20), (29), and (34). 
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The flows represented by equation (34) are actually a special case 
of the exact solutions of reference 15 and will not be discussed further. 
As a result of these analyses, table I has been prepared, which summarizes 
the four new cases obtained here of mainstream flows over a flat or 
slightly curved surface for which the boundary- layer flows have similarity 
with respect to polar coordinates. 

The cases are as follows: 

Case I: U = ar n e m+1 

W = br n 0 m , b ^ 0 
Case II: U = ar n 0 m+1 

W = 0 

W = br n 0 m , b ^ 0 
Case III: U = af(r)0 m_1 

W = bf(r)0 m , b ^ 0 
Case IV: U = af(r)0 m_1 

W = 0 

W = bf(r)0 m , b ^ 0 

As described earlier, secondary flows exist even though the radial 
component of mainstream flow vanishes (W = 0). For such cases, a function 
W = W(r ,0) ^ 0 is defined, and the boundary- layer radial component of 
flow is expressed as 


w = WG * ( T] ) 


(4c) 


The Mainstream 

When W ^ 0, the mainstreams are spiral flows. For W = 0, circular 
mainstream flows are obtained. 

U, W, W. - In regions where the thin-boundary- layer theory is ap- 
plicable, the mainstream is very nearly parallel to the surface; U and 
W are functions of r and 0 only. 

The analysis for the exact solutions for boundary- layer flows having 
similarity with respect to stationary polar coordinates (ref. 15) showed 
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that only one form of U and W (or w) is possible, that is 

U = ar n e me 

fl (36) 

W (or W) = br n e m 

The present analysis results in obtaining a much wider variety of flow 
solutions, represented by equations (20) and (29). In particular, there 
are no restrictions whatsoever on the functional relation of U and W 
(or W) to r in cases III and IV (eq. (29)). An example might be 
chosen of circular main-flow streamlines where the inlet velocity U at 
the leading edge 0=0 varies in a sinusoidal fashion. The secondary- 
flow overturning would then cause the boundary layer to pass through 
regions where the main-flow velocity is alternately increasing and de- 
creasing. As will be seen later in the discussion of the boundary- layer 
thickness and g(r,0), the boundary layer will correspondingly become 
thinner or thicker. 


The small- angle solutions were obtained by assuming that second and 
higher power terms of 0 could be neglected. In all these cases, U and 
W (or W) are obtained as products of powers of r and 0. In cases I 
and II W/U (or W/u) is proportional to l/d. In cases III and IV, w/U 
(or w/u) is proportional to 0. 

Projection of main-flow streamline on surface . - The equation for the 
projection of the main- flow 1 streamline on the surface for W / 0 may be 
obtained by integrating 


W _ dr 
U ~ r d0 


(37) 


Whenever W / 0 and U / 0, spiral mainstream flows result (cases I and 
III). For W = 0, circular main- flow streamlines result (cases II and IV). 

Slope of projected streamlines. - The slope of the projected stream- 
line with respect to 6 = 0 (the tangent of the angle between the tangent 
to the projected streamline curve at a point and the line 0=0) may be 
obtained from 


slope = 


dr 

^ tan 9 + r 

U.C7 

- r tan 0 

UC7 


(38) 


It is found by substitution into equation (38) that the slope is inde- 
pendent of radial position r for all cases. 
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Irrotat ionality . - For mainstream flows considered here and in re- 
gions of thin boundary layers, as required for this analysis, only the 
component of vorticity normal to the surface 


1 dW _ dU _ U 
r c(0 3r r 


(39a) 


can be much different from zero (ref. 14). The values of the constants 
specified under the listing "Irrotationality" (table i) were obtained in 
each case from 


1 dW dU _ U 
r ^0 ^r r 

These values set the conditions for nearly irrotational mainstream flows. 



The Boundary Layer 

As discussed in references 14 and 15, the physical interpretation of 
the boundary- layer behavior that the mathematical representations describe 
is best found by examining the behavior of t) and in particular g(r , 9 ). 

The boundary- layer thickness on the surface at a point r,0 is in- 
versely proportional to g(r,0) at the point. In order for the theoreti- 
cal boundary layer to have a beginning at a leading edge with zero thick- 
ness, as in a real fluid, there should be a line along the surface for 
which g(r,0) is infinite while the velocities remain finite. In the 
exact solutions presented in reference 15, this occurs in the finite part 
of the plane only at the point r = 0 for values of n < 1. For n > 1 
the boundary layer in reference 15 may be considered to have a "beginning" 
only at r = ». However, the mainstream velocities there take on 
" infinite" values . 

In the present report (flows for small angle 0), for cases I and II 
when -1 5 m < 0, g(r ,0) and hence rj becomes infinite along the line 
0=0, while U is finite as required. In case I, however, W (and w) 
are unbounded, so a proper leading edge does not exist there. In case 
II, W = 0 but W is unbounded at 0=0. In cases III and IV, for 
1 < m <2, there is a properly defined leading edge at 0=0, where 
g(r,0) takes on infinite values and where U, u, W, and w remain 
finite. For 1 < m < 2, in cases III and IV, U = W = 0 at the leading 
edge . For m = 1, U is independent of 0 and may be different from 
zero at the leading edge . 
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In cases I and II for m > 0, and in cases III and IV for m > 2, 
along the line 8=0, the mainstream and the boundary- layer velocity 
components are all zero. Even though the boundary- layer velocities match 
the mainstream in these situations, a proper leading edge does not exist 
there because g(r,0) equals zero, corresponding to an "infinitely" thick 
boundary layer. As in reference 14, such accelerated- flow cases may be 
considered appropriately by confining the discussion to regions where the 
requirement of thin boundary layers is satisfied. 


The Ordinary Differential Equations 

The actual numerical solutions of the ordinary differential equa- 
tions are beyond the scope of the present investigation. The literature 
contains examples of numerical solutions for particular values of the 
constants. Some of these examples are noted in the listing "Comments and 
References" associated with each case in table I. 

The present analysis simply derives the ordinary differential equa- 
tions that can be obtained with the underlying assumptions. In any 
particular case of interest for which the equations are appropriate, the 
existence of the numerical solution and its computation must be obtained 
individually. Nevertheless, some general remarks (in part repeating 
material from ref. 14 here for convenience) can be made here (as in ref. 

II) concerning the numerical solutions. 

Separation of F and G . - Under certain choices of the free con- 
stants involved, the functions F and G are separable; that is, one 
equation of the pair of ordinary differential equations will contain 
terms in only one of these functions and its derivatives. Numerical 
solutions are much more readily obtained in such cases than when the 
functions are not separated. 

An example is provided when a = 0 (case i), so that equation (22) 
then contains terms only in G and its derivatives. In all cases where 
the functions can be separated, the equation is thereby reduced to a 
Falkner-Skan type equation. The complete solutions to Falkner-Skan equa- 
tions have been obtained in references 4 and 5. Thus, equation (30) (case 

III) and equation (32) (case IV), in which the functions are already 
separated, are all Falkner-Skan equations. 

Although it is not apparent from the equations in the table alone, 
when a = 0, then u = 0 and equation (la) and therefore equations (2l), 
(23), (30), and (32) disappear. In cases I and III (W ^ 0), when a = 0, 
the flows are straight two-dimensional flows along radial lines out from 
a stagnation point. In cases II and IV (W = 0), a=0 is the trivial 
case of no mainstream flow. 
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In case III, when a = b, F = G. If, in addition, m = 0 and 
c = -a, equation (30) "becomes equation (9.12) of reference 16: 


As pointed out in reference 16, this is one of the few cases when the 
boundary- layer equation can be solved in closed form. The solution is 


then by differentiation of equation (43) and substitution, equation (42) 


Pij^^girity in u or w. — As discussed in reference 14 and applied 
in reference 13, an extension of the solutions beyond strict similarity 
of the velocity component can sometimes be made by addition of solutions 
where the boundary- layer equations are linear in u or w. Apparently 
such extensions are not possible for the boundary- layer flows investi- 
gated here, because equation (l) is always nonlinear in u and in w 
except for the typical case of no mainstream flow. Nevertheless, there 
is the complete freedom from specification of the form of f(r) in cases 
III and IV. Accordingly, if 


F'" - F‘ 2 +1 = 0 


(40) 



(41) 


When c = a = b and m = 0, equation (30) becomes 

pm + f* 2 - 1=0 


(42) 


Letting 


= F(n) 


(43) 


becomes 


jr»* - Sf' 2 - + 1 = 0 


(44) 


and the solution thereby is seen to be 



(45) 


U 1 = a 1 f 1 (r)0 m_1 


U 2 = a 2 f 2 (r)e m - 1 


and 
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are possible solutions for mainstream flows, then 


U s U-l + U 2 


fl f l( r ) + a 2 f 2 (r)j e ™' 1 s af(r )0 m_1 


likewise represents a possible mainstream flow, although solutions for 
u and w cannot be superimposed. This fact (the reasoning is the same 
for W as for u) is evident also from the ordinary differential equa- 
tions for cases III (eqs. (30) and (3l)) and IV (eqs. (32) and (33)), 
which are independent of the form of f(r). 


An experimental investigation was made to determine whether the 
theory provides a reasonable approximation of the limiting flow deflec- 
tions for the particular case of circular flow over a flat plate. Of the 
cases presented, this case is the one most likely to be encountered in 
actual practice. 

Theoretical prediction of limiting deflection . - The case that will 
be investigated is case TV with m = 1. Under this assumption equations 
(32) and (33) become, respectively, 


where the following relation between the various constants has been 
chosen: 


Equation (46) is the Blasius equation, and values for F are tabulated 
in reference 13. Reference 12 shows that the function G 1 in equation 
(47) is expressible in the form 


Comparison with Experiment 



(46) 



(47) 


c = a = -b 


G« = P(n) - F' (t| ) 

where P(tj ) is the solution of the equation 


(48) 



( 49 ) 
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with the boundary conditions 

lim P(tj) = 1, P(0) = 0 


The solution of equation (49), however, is presented in reference 13, and 
the values of P(tj), (P^t)) in ref. 13) are presented in tables. 

Now the angle of flow deflection of the boundary layer at the plate 
surface is determined by 

X = arc tan f lim — ) 

Uo u / 

= arc tan ( lim \ 

v» / 

As G'(0)/F*(0) is an indeterminate form, application of L' Hospital's 
rule gives 

T = arc tan ( lim |rr ) 

\1 4 0 / 


= arc tan I lim 
V! 0 


pi _ pn\ 


F" 


= arc tan 


lim ff£ - 1 
T) ■+ 0 \ 


The value of P'(0)/F"(0), however, is determined in reference 12 to be 
4.270. Hence, 


T = arc tan 3.270 


(50) 


From equation (50), the equation for the limiting deflection line on 
the plate surface can he found as the solution of the differential 
equation: 


- ~ = tan r = 3.270 
r d0 


(51) 


Equation (5l) has the solution 


3.270 9 


r = (const. ) e' 


(52) 
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Experimental determination of limiting deflection . - The experimental 
determination of limiting flow deflection was made by means of smoke flow 
visualization (apparatus and procedures described in ref. l) in a Lucite 
two-dimensional circular channel of rectangular cross section. Tests were 
conducted on a plate parallel to the base of the channel and fastened to 
the channel walls at a distance of approximately one-third the channel 
height. With a maximum Reynolds number of about 6><10^, the flow for these 
tests was well within the laminar range. 

Figure 1, a view through the Lucite top of the channel, shows smoke 
introduced into the mainstream at approximately a midchannel position. 

On the test plate shown in the photograph, the circular white lines cor- 
respond to the theoretical mainstream flows . The straight lines orthogo- 
nal to the circular lines represent 10° increments in 9. The dotted 
lines appearing on the plate are the theoretical limiting deflection lines 
based on equation (52). The main- flow streamlines as depicted by the 
smoke trace were found to follow closely the theoretical main- flow stream- 
lines well beyond the region of interest near the leading edge of the 
plate, although this is somewhat obscured by parallax in figure 1. Near 
the exit of the channel, some deviation of the smoke trace from the cir- 
cular arcs occurred because of secondary- flow accumulations. 

Figures 2(a) and (b) show limiting flow deflection determined by 
introducing smoke directly on the plate surface. The photographs show 
that the theory predicts limiting deflection very well in the range 
0=0^ 10°. Beyond 10° the theory predicts a greater overturning of 
flow than that indicated by the experiment. These results give an order 
of magnitude to the range of values of 9 where the small- angle approxi- 
mation appears to be reasonable. 

Prediction of limiting deflection based on translate flow . - Refer- 
ence 12 presents an analysis of three-dimensional boundary- layer flows 
when the main- flow streamlines are translates (i.e., all streamlines are 
obtained from a single streamline by propagation of the streamline par- 
allel to the plate leading edge). As a matter of interest, this theory 
was also used to predict limiting deflection for the present case. 

In effect, a comparison can be made between the results predicted by 
two distinct kinds of approximations. In the small- central- angle varia- 
tion method the approximation is made in the solution of the boundary- 
layer equations while the description of the main flow is exact. In the 
streamline- translate method the approximation is made in the representa- 
tion of the main-flow streamlines while the solution to the boundary- layer 
equations obtained is exact. 

For the streamline- translate method, a mean value between the radius 
of the inner wall of the channel and the radius corresponding to the 
starting point of the outer limiting line shown plotted in the photographs 
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was chosen as the defining radius of a typical circular-arc streamline in 
translate flow. Since circular-arc translate flow is analyzed and the 
limiting flow deflection is determined in reference 12 (see fig. 5, ref. 

12 ), the results could be applied directly in the present investigation. 
The basic assumption underlying such an application is that circular-arc 
translate flow of the type described is a reasonable approximation of the 
concentric-circular-arc streamline flow that actually exists in the 
channel. 

The results of the analysis are presented in figure 2(c). A 30° sec- 
tion of the channel is sketched, and the theoretical limiting deflection 
lines for both the translate- flow analysis and the small-angle analysis 
are indicated. The circles appearing on the sketch represent points on 
the smoke flow pattern as measured from the photographs. It can be seen 
from figure 2(c) that the translate- flow analysis is less accurate than 
the small- angle analysis up to about 10° •, but, from 10° on, the agreement 
between theory and experiment improves for the translate analysis. After 
20° it appears that the predicted limiting line and the actual limiting 
line are very nearly parallel up to the point where wall interference 
causes the boundary layer to deflect in a circumferential direction. 


CONCLUDING REMARKS 

Exact solutions were obtained in reference 15 describing the main- 
stream flows over a flat or nearly flat surface for which the thin laminar- 
boundary- layer flows have similarity with respect to stationary polar co- 
ordinates. The solutions thus obtained were of the form 

U = ar n e m0 

W = br n e m0 

By suitable choice of the constants a, b, n, and m, the main flows may 
be stagnation- point, spiral, or circular flows. The boundary layers for 
these exact solutions have no properly defined leading edge in the finite 
part of the plane, resulting in some awkwardness in relating these theo- 
retical flows to real physical flows. 

The present analysis is restricted to regions of small central 
angle 0. The appropriateness of the restriction in physical situations 
has been established experimentally (refs. 1 to 3 and 13). The mainstream 
flows obtained herein are 

TI nun+1 

U = ar 6 

W (or W) = br n 0 m 


/ 


(cases I and II ) 
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U = af(r)0 m_1 
W (or W) = bf(r)0 m 


(cases III and IV) 


In cases III and IV for 1 < m < 2, a properly defined leading edge can 
be obtained along the line 0=0 where the boundary layer has zero . 
thickness and the mainstream velocity components do not become infinite. 
Although different velocity distributions are obtained here, the projected 
main-flow streamline configurations possible are the same as those of ref- 
erence 15; that is, (l) stagnation flows along radial lines from a stag- 
nation point, (2) spiral flows out from (or in toward) a central point, or 
(3) circular flows. 


Actual numerical solutions of the transformed boundary- layer equa- 
tions are not attempted here. Particular examples are noted, however, 
for which solutions have been obtained elsewhere. 

An experimental comparison of limiting overturning at the wall under 
circular mainstream flow using smoke flow- visualization techniques showed 
good agreement for the first 10° of the mainstream turning. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, September 26, 1956 
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TABLE I. - SIMILARITY SOLUTIONS IN STATIONARY POLAR COORDINATES FOR SMALL ANGLE VARIATION 


(a) Case I 


u 

n m+l 
ar 9 

w 

, n _m , , ^ 

br 0 , b / 0 

V 

/cr n 1 0 m N / cU \ /cW \ 

V J y\var&J - y^vbrj » c / 0 

Ordinary 

differential 

equations 

(21) b (n + 1)(F'G' - 1) - b -^ n ^ 5 ) GF" + a(m + l) [(F 1 ) 2 - l] - a ^ m g +2 ) FF" - cF'" = 0 

(22) am(F'G' - l) + bn [jG 1 ) 2 - l] - — GG" - a ^ m g 2 ^ FG" cG'" - 0 

Boundary conditions 

F 1 (0) = G'(0) = F(0) = G(0) = 0, lim F ' (7?) = lim G'(7?) = 1 

77 -► 00 77 °° 

Projection of main- 
stream on surface 

r = C0 b//a spiral flow streamlines (a / 0) 

Slope of projected 
streamline with 
respect to 0=0 

— tan 0+0 
a 

- 0 tan 0 
a 

Irrotationality 

m = a(n + l) = 0 

Linearity in 
u and w 

[lb | possible for this case 

Separation of 
F and G 

(21) Not possible for this case 

(22) a = 0 

Comments and 
references 

a = 0, stagnation flow: eq. (la) and eq. (21) vanish and eq. (22) becomes a 

Falkner-Skan equation with solution completely known, refs. 4 and 5. 

Ref. 8: a = c = n = -1, m = b , b = -c from ref. 8. 
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TABLE I. - Continued. SIMILARITY SOLUTIONS IN STATIONARY POLAR COORDINATES 


w 

o 


FOR SMALL ANGLE VARIATION 
(b) Case II 


u 

ar n e m+1 

w 

0 

w 

br n 0 m , b / 0 

17 


Ordinary 

differential 

equations 

(23) b (n + 1) F'G' - GF " + a(m+l)[(F') 2 - 1] - a(m + 2 ) FF" - cF"’ = 0 

(24) am(F'O') + bn(G') 2 - b < n £ 5 >- GO" - FG” - CG>" - 0 

Boundary conditions 

p* (0) = G ' ( 0 ) - F(0) = G(0) = 0, lim F • (77) = 1, lim G'(7?) = 0 

7?-* o° 7? -► ® 

Projection of main- 
stream on surface 

r = C, circular-flow streamlines 

Slope of projected 
streamline with 
respect to 0=0° 

- Cotangent 0 

Irrotationality 

(n + 1) =0 

Linearity in 
u and w 

[lb| Not P osslble for thls case 

Separation of 
F and G 

[if] Not P ossil:>le f° r this case 

Comments and 
references 

a = 0, no flow 
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TABLE I. - Continued. SIMILARITY SOLUTIONS IN STATIONARY POLAR COORDINATES 

FOR SMALL ANGLE VARIATION 


(c) Case III 


u 

af(r)e m_1 

w 

bf(r)6 m , b / 0 

7? 

y (klz^) =y (^) = y ( v - 2 ) .c/0 

Ordinary- 

differential 

equation 

(30) a(m - 1) |jF-) 2 - l] - FF" - cF"' = 0 

(31) - -^[(F 1 ) 2 - l]+ am ( F 1 G ' - l) - FG" - cG'" = 0 

Boundary conditions 

F ' (0) = G ' (0 ) = F(0) = G (0) = 0, lim F ' ( 77 ) = lim G ' ( 7? ) = 1 

Projection of main- 
stream on surface 

r = Ce b ^ 2 / 2a , spiral flow (a / 0) 

Slope of projected 
streamline with 
respect to 0=0 

0 tan 0 + ^ 
0 - ^ tan 0 

Irrotationality 

£ bf (r) - af'(r) - | f(r) = 0 

Linearity in 
u and w 

Not possible for this case 

Separation of 
F and G 

(30) Separated 

(31) a =* m = 0 (see comments and references) 

Comments and 
references 

Eq. (30) is a Falkner-Skan equation, completely solved 
in refs. 4 and 5. 

If a = 0, eq. (la) and eq. (24) vanish, stagnation flow, 
boundary conditions not achievable. 

If a » b, F » G: if m « 1, F = Blasius F 

Ref. 16, eq. (9:12): a = b = -c < 0, m = 0. 
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TABLE I. - Concluded. SIMILARITY SOLUTIONS IN STATIONARY POLAR COORDINATES 

FOR SMALL ANGLE VARIATION 
(d) Case IV 


Ol 

CO 


u 

af(r)e m_1 

w 

0 

w 

bf(r)e m , b / 0 

T) 

• / o 

Ordinary 

differential 

equation 

(32) a(m - l) |[F' ) 2 - l] - FF" - cF'" = 0 

(33) - ^ [(F 1 ) 2 - l] + am ( F ' G ' ) - -^r FG" - cG'" » 0 

Boundary conditions 

F' (0) = G' (0) = F(0)= G(0) = 0, lim F 1 (77 ) = 1, lim G'(l))=0 

r?-*~ 

Projection of main- 
stream on surface 

r = C, circular-flow streamlines 

Slope of projected 
streamline with 
respect to 0=0 

- Cotangent 0 

Irrotationality 

f'(r) + i f(r) = 0 

Linearity in 
u and w 

[lb] P° ss itl e for this case 

Separation of 
F and G 

(32) Separated 

(33) a = 0 (See comments and references) 

Comments and 
references 

Eq. (32) is a Falkner-Skan equation completely solved in 
refs. 4 and 5. 

If a = b, F ^ G (see boundary conditions) 

If m = 1, F = Blasius F 
Ref. 9: c = a = b, m = 1. 

If a = 0, no flow. 
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(a) Boundary- layer streamline near pressure surface at inlet. 
Figure 2. - Limiting flow deflection in circular channel. 
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(h) Boundary- layer streamline near midchannel at inlet. 

Figure 2. - Continued. Limiting flow deflection in circular channel. 
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(c) Comparison of translate solution with small-angle solution and smoke pattern. 
Figure 2. - Concluded. Limiting flow deflection in circular channel. 
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